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1. Introduction 

Calculating nonperturbative properties of supersymmetric theories on the lattice encounters 
various difficulties related to the fact that the discretisation of space-time explicitly breaks super- 
symmetry and violates Leibniz' rule. Moreover, the vanishing of the Witten index in the context 
of spontaneous supersymmetry breaking leads to a fermion sign problem which makes straightfor- 
ward numerical simulations impossible. While the restoration of supersymmetry can sometimes be 
achieved in the continuum limit of the lattice theory, e.g. by fine tuning or by constructing Q-exact 
discretisations [1], a solution to the sign problem is not easy to find. A possible way out has been 
proposed in [2, 3, 4]. It is based on the fermion loop formulation which can be simulated without 
critical slowing down even when a massless goldstino mode is present. 

In this work, we apply the fermion loop formulation to JV = 2 supersymmetric quantum me- 
chanics for superpotentials yielding broken or unbroken supersymmetry. Using transfer matrix 
techniques we are able to obtain exact results for partition functions and various observables at fi- 
nite lattice spacing. We investigate how the supersymmetric spectrum is recovered in the continuum 
limit and how the goldstino mode emerges in the case of broken supersymmetry. In these proceed- 
ings we confine ourselves to the presentation of results obtained using a Wilson type discretisation 
together with the appropriate fine tuning of counterterms, although results using a Q-exact discreti- 
sation have been derived as well. 



2. Supersymmetric quantum mechanics on the lattice 

The continuum action of ,jY = 2 supersymmetric quantum mechanics can be written as 

1 (m) ^ ''+\p'm)) 2 +w(t)[j t +p"(m) v(t) 



dt 



2 \ dt 



(2.1) 



with one real bosonic coordinate 0, two anticommuting fermionic coordinates y and i/f, and a 
generic superpotential P((p). The derivative of the superpotential P'(<p) is taken with respect to 0, 
P'((j)) = dP g^ . For periodic boundary conditions (PBC) the action is invariant under two super- 
symmetry transformations 8^: 

Slip = 1/A£, #20 = 



8 lV = 0, 52 ¥ =(f-P')e, 
8 1 y=(f+P>)e, 521/7 = 0, 



with two Grassmann valued parameters e and e. Note that for supersymmetric quantum mechanics 
it is the form of the superpotential P(<j>) which determines the supersymmetry breaking pattern. If 
the highest power of P(</>) is even (odd), supersymmetry is unbroken (broken). A main feature of 
supersymmetry is the degeneracy between the energy levels in the bosonic and the fermionic sector. 
For unbroken supersymmetry, however, there is one single unpaired energy level at zero energy, i.e., 
a unique ground state, either in the bosonic or in the fermionic sector. This is in contrast to the case 
of broken supersymmetry, where the lowest energy levels in both sectors are degenerate and lifted 
above zero. In addition, there is a zero energy goldstino mode which mediates between the two 
degenerate ground states. 
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The supersymmetry breaking pattern can also be partly infered from the Witten index. It is 
formally defined as 

W = lim Tr [(-l) F exp(-/3#)] , 

where F denotes the fermion number operator and H is the Hamiltonian of the system. Essentially, 
W counts the difference between the number of bosonic and fermionic zero energy states and its 
vanishing provides a neccessary but not sufficient condition for supersymmetry breaking. W can 
also be written more explicitely as 

W = lim [Tr b exp(-pH) -Tr/exp(-j3#)] = lim [Zq-Z{\ = lim Z PBC , (2.2) 

where Tr&j denote the traces over the bosonic and fermionic states. Zoj are the partition functions 
in the F = 0, 1 sectors and Zpbc is the one with periodic boundary conditions. In the language of 
field theory, the latter can be calculated via 

Z PBC = J @<S>9\jr3iyexv{-S) = J 00detD(0)exp(-fy). 

In the last step, the fermions have been integrated out yielding the fermion matrix determinant 
detD and the bosonic part of the action S$. In this representation, the origin of a fermion sign 
problem becomes evident when supersymmetry is broken: a vanishing Witten index requires the 
determinant detD to be indefinite. 

2.1 Lattice formulation 

For the construction of a lattice version of the model, we follow Golterman and Petcher [5] and 
employ the same lattice derivative for the bosons as for the fermions. To avoid fermion doublers, 
we use the Wilson lattice derivative with Wilson parameter r = 1 . In one dimension this simplifies 
to the backward derivative (A~f) x = f x — f x -\ and the discretised action explicitly reads 



~(P'(^) 2 +2^) - + (1 +P"(4>x))¥xYx ~ Wx¥x-l 



(2.3) 



Due to radiative corrections the lattice theory is, however, not guaranteed to yield a supersymmet- 
ric theory in the continuum limit. The corrections can be accounted for either by adding a suitable 
counterterm \ Y*P" to tne action [5,6], which restores the supersymmetries in the continuum limit, 
or by adding the surface term £P'(A~0) [7, 8, 9] resulting in a Q-exact action. The latter construc- 
tion preserves a particular combination of the supersymmetries Sip exactly even at finite lattice 
spacing and hence guarantees the correct continuum limit without any fine tuning. 

To circumvent the sign problem discussed above, we make use of the fermion loop formulation 
[2, 3, 4]. The basic idea here is to exactly rewrite the exponential of the fermion degrees of freedom 
as a power series to all orders. Upon integration of the fermion fields, the nilpotency of the Grass- 
man variables yields a constraint on the oriented fermionic bond occupation numbers n x = 0, 1 
between the sites x and x— 1 related to the fermion hopping term \j/ x y x ^i, and on the monomer 
occupation numbers m{ = 0, 1 stemming from the term (1 + P" (^ x ))\j7 x Yx- The constraint is given 
by 

m {+~{n x *+ nf x+i) = 1 Vx > 
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and allows only two fermion configurations: {m x = i,n x = 0, Vx} with fermion number F = 0, 
and {m{ = 0,n£ = 1, Vx} with fermion number F = 1. For PBC the latter receives an additional 
minus sign relative to the former due to the fermion loop. As a consequence, the partition function 
naturally decomposes into a bosonic and fermionic contribution Zo and Z\, in accordance with 
eq.(2.2). It is this decomposition which eventually allows to take care of the fermion sign problem. 

In addition to the fermion bonds and monomers we also introduce non-oriented bonds for the 
bosonic degrees of freedom, with the corresponding bosonic bond occupation numbers n b € N° 
[4]. Depending on the symmetries of the action, these bosonic bond configurations may obey 
certain constraints. By summing over these constrained configurations {n b } we obtain the locally 
factorised partition functions with fixed fermion number F = 0,1, 

where the local weights Qp are defined as 

Q F {N) = Jd<$> w e4(^)W) (1+jP - (0)) i-F 

with the bosonic site occupation number N x = n b + n b +l . The <2-exact discretisation requires addi- 
tional types of bosonic bonds, but still leads to a locally factorised partition function. 

2.2 Transfer matrix 

The dimensionality of the system allows a further reformulation in terms of a transfer matrix 
between states defined on the dual lattice. Each state is characterised by the fermion bond occupa- 
tion number number nf and the boson bond occupation number n h , i.e. \n? ,n b ). Since the fermion 
number is conserved the transfer matrix has a block structure consisting of the two matrices T Z~ h ' 

O m tt n B 

which take the system from state \F,n b ) to \F,m b ). To be specific, the transfer matrix elements are 
given by 

T F b b = ^=^=Q F (m b + n b ). 

m°,n / uT v ' 

Vm D \ Vn°\ 

In order to keep the size of the matrices finite, we introduce a cutoff on the maximal bosonic bond 
occupation number. Keeping it of the order <^(10 2 ) turns out to be sufficient to render all results 
independent of the cutoff. 

In terms of these transfer matrices, the partition function for a system with L t lattice sites is 
calculated in each sector F according to 

Z F =Tr[(r F ) L ']. 

These partition functions can then be combined to Zpbc = Zq—Z\ and Z u pbc = Z$+Z\ for PBC and 
antiperiodic boundary conditions (aPBC), respectively. The construction via the transfer matrices 
allows the straightforward calculation of various observables, such as correlation functions, Ward 
identities and mass gaps. The latter are directly associated with the eigenvalues of the transfer 
matrices. If we denote the eigenvalues of T F by Xq > Af > . . . , the k-th bosonic mass gap in the 
sector F can be calculated as 

m F b k = -L r \og(Xi/Xl), k=l,2,..., 
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Figure 1: The Witten index W = ZpBc/Z a pBC versus the lattice spacing for broken (left plot) and unbroken 
supersymmetry (right plot) at various values of the inverse temperature mL. The continuum limit corresponds 
to am — > 0. 



whereas the k-th fermionic energy gap is given by 

mf = -L r log(Xl- F /tf), k = 0,l,.... 

The transfer matrix approach can, of course, be generalised straightforwardly to any kind of dis- 
cretisation of the action eq.(2.1), in particular also to the Q-exact discretisation. 



3. Results 

We now present the results for the action (2.3) with the counterterm using the techniques 
introduced above. For our calculations, we use the superpotential P u ((j)) = ^m(j) 2 + \g<j> 4 as an 
example with unbroken supersymmetry and Pb(^>) = — ^0 + ^A0 3 as an example for which the 
supersymmetry is broken. The calculations are performed at coupling strengths g/m = 1.0 and 
X jm > l 1 = 1.0, respectively, thus we are clearly in a regime where perturbation theory is not appli- 
cable. For a system with aPBC for the fermion the temporal extent of the lattice is inversely related 
to the temperature T of the system, such that mL — > oo corresponds to the zero temperature limit. 
Finally, the continuum limit is reached by taking L t — > oo. 

3.1 Witten index 

The Witten index is determined by the quantity Zpbc/ZuPBC = (Zo —Z\)j (Zo +Z\ ). It measures 
the relative weight between the bosonic and fermionic sectors Zo and Z\ , respectively. In the system 
with broken supersymmetry both ground states are equally favourable, yielding ZpBc/Z a pBC = 
in the zero temperature limit. Of course, the degeneracy between the two ground states is broken 
at finite lattice spacing, so one expects a Witten index W = ±1 in the limit T — > at fixed a. It 
turns out that for our choice of parameters, the fermionic ground state has a slightly lower energy 
at finite a leading to W = — 1 in the T — > limit, cf. left plot in figure 1 . This is true for any finite 
a, so the order of the limits lim^o lim a ^o is crucial to obtain W = 0. Note also that for T S> 1 the 
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Figure 2: Bosonic and fermionic mass gaps mi, and nif versus the lattice spacing a, all expressed in units of 
the bare mass m. Results for broken and unbroken supersymmetry are displayed in the left and right panel, 
respectively. Note that for broken supersymmetry the zero energy Goldstino mode emerges in the continuum 
limit. 



Witten index tends to zero at any finite a, in accordance with the counting of the states in eq.(2.2) 
at finite temperature. 

In the situation with unbroken supersymmetry the system is forced to occupy the single unique 
ground state in the zero temperature limit, yielding Zpbc/Z c ,pbc = +1 or — 1. Note that for our 
specific choice of parameters the ground state is bosonic, hence W = +1. However, while one 
finds that the index is pushed away from for T S> 1 as before, in the limit T — > it will always go 
to 1 at any finite a, cf. right plot in figure 1. So it turns out that for unbroken symmetry, the order 
of the two limits lim^o and lim u ^o is not relevant. 



3.2 Mass gaps 

It is also interesting to study how the energy or mass gaps approach the continuum limit. In 
figure 2 we show the results for the lowest few masses as a function of the lattice spacing a, every- 
thing expressed in units of the bare mass m, for broken (left plot) and for unbroken supersymmetry 
(right plot). Since we extract the mass gaps from the eigenvalues of the transfer matrices, the re- 
sults are obtained directly in the limit T — > 0. As a consequence, for broken supersymmetry, where 
there are two degenerate ground states in the continuum, it makes sense to calculate bosonic exci- 
tations m/, both in the F = and F = 1 sector. The plots illustrate nicely how the supersymmetry 
in the spectrum, i.e. the degeneracies between the bosonic and fermionic excitations, are restored 
in the continuum limit. Furthermore, when the supersymmetry is broken one expects a zero en- 
ergy fermionic excitation, the goldstino mode, which is responsible for the fact that Zpbc = 0. 
From the plot it becomes clear how the lattice acts as a regulator for the goldstino mode and, as 
a consequence, also for the vanishing Witten index W, hence allowing to give meaning to (finite) 
observables even in the system with PBC. Finally, we make the observation that the leading lat- 
tice artefacts of the spectral mass gaps are all &(a) except for in the F = sector when the 
supersymmetry is broken. In that case they are &(a 2 ). 
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4. Conclusions 

We have presented exact results for JV = 2 supersymmetric quantum mechanics on the lattice 
using the fermion loop formulation and corresponding transfer matrices. With these techniques we 
are able to study in detail how the supersymmetric spectrum is recovered in the continuum limit 
and how the Witten index is regularised on the lattice. 

In the loop formulation the partition function naturally separates into bosonic and fermionic 
contributions and this is crucial for containing the fermion sign problem in supersymmetric sys- 
tems with broken supersymmetry. The transitions between the bosonic and fermionic sectors are 
controlled by the (would-be) goldstino mode which becomes massless only in the continuum limit. 
Since massless fermion modes can be efficiently simulated with the fermion loop algorithm pro- 
posed in [2, 3] our approach provides a way to circumvent the sign problem. Indeed, results from 
Monte Carlo simulations of jY = 2 supersymmetric quantum mechanics have already been pre- 
sented in [4] and in this work we have provided the corresponding exact results using transfer 
matrices. 

It is also interesting to apply our approach in higher dimensions where it allows to investigate 
the spontaneous breaking of supersymmetry nonperturbatively and from first principles. In partic- 
ular, the approach can be applied to supersymmetric Wess-Zumino models [4] in d = 2 dimensions 
and first results from simulations of the jY = 1 model including one Majorana fermion and one 
scalar field have been presented at this conference [10]. 
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